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The onset of generalized synchronization of chaos in directionaUy-coupled systems corresponds to 
the formation of a continuous mapping which enables one to persistently define the state of the 
response system from the trajectory of the drive system. The recently developed theory of gener- 
alized synchronization of chaos deals only with the case where this synchronization mapping is a 
single-valued function. In this paper, we explore generalized synchronization in a regime where the 
synchronization mapping can become a multi-valued function. Specifically, we study the properties 
of the multi- valued mapping which occurs between the drive and response systems when the systems 
are synchronized with a frequency ratio other than one-to-one, and address the issues of the exis- 
tence and continuity of such mappings. The basic theoretical framework underlying the considered 
synchronization regimes is then developed. 
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I. INTRODUCTION 

The synchronization of oscillations is one of the most 
interesting nonlinear phenomena and is an inherent part 
of many processes studied in a wide range of natural 
systems, including such diverse areas as neuro-biological 
networks and solitary systems. As a result, the corre- 
sponding theory is extensively utilized in many practical 
applications [Q-^ ■ The dynamical theory of the synchro- 
nization of periodic oscillations, which relates the onset 
of synchronization to the birth of a stable limit cycle due 
to the bifurcation of the motion on a two-dimensional er- 
godic torus, is due to Van der Pol [|j . Since the discovery 
of chaotic behavior in nonlinear oscillators and the abil- 
ity of chaotic oscillators to synchronize, the framework of 



the dynamical theory of synchronization has been signif- 
icantly modified. The theory now encompasses the ma- 
jor properties inherent in the synchronization that deals 
with limiting sets (called chaotic attractors) which are 
dynamically much richer than isolated limit cycles |^]. 

Different notions of synchronization occurring in 
chaotic oscillators have been introduced in order to ex- 
plore the qualitative dynamical changes caused by the 
onset of chaos synchronization in numerous experimen- 
tal studies. Such notions include the cases of identical 
synchronization, where identical chaotic oscillations are 
usually studied in coupled systems with identical indi- 
vidual dynamics [^j^, (7enerafeed synchronization, which 
extends the notion of identical synchronization to cases 
of directionaUy-coupled systems with non-identical in- 
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dividual dynamics [^|J^-[ll|], and phase synchronization 
which usually deals with the phase locking of the main 
frequencies in the spectrum of the chaotic systems while 
the chaotic components of the signals remain indepen- 
dent [ p^yi^ ]. Studies of the phenomena of generalized 
and phase synchronization are largely motivated by the 
need for development of a theoretical framework that 
can be used to achieve a better understanding of syn- 
chronization in experiments with neuro-biological sys- 
tems where complex, chaotic behavior of neurons is very 
typical gjlig. 

This paper addresses issues of generalized synchroniza- 
tion which arise in studies of chaos synchronization with 
a frequency ratio other than one-to-one. As it is defined 
in 1^ , the onset of generalized synchronization in direc- 
tionally coupled chaotic systems relates to the formation 
of a continuous mapping that transforms the trajectory 
on the attractor of the drive system into the trajectory 
of the response system. For systems with invertible dy- 
namics, this is equivalent to the formation of a contin- 
uous mapping which links the current states of the sys- 
tems once they settle down on the synchronous attractor. 
Since the introduction of the concept of generalized chaos 
synchronization, significant progress has been made in 
understanding the relationship between the properties 
of these synchronization mappings and the spectrum of 
Lyapunov exponents which characterize the synchronous 
chaotic attractor. A number of papers have analyzed con- 
ditions which guarantee the differentiability of the syn- 
chronization mapping, thus indicating the formation of a 
differentiable invariant manifold (in the joint phase space 
of the coupled systems) which contains the synchronized 
chaotic attractor (see, for example [|r^,^). Simply put, 
this differentiability occurs when the rate of contraction 
in the direction transverse to the manifold is larger than 
the rate of contraction experienced by trajectories of the 
chaotic attractor in the direction tangent to the manifold. 
Such cases are distinguished as belonging to a subclass 
of generalized synchronization which has been named dif- 
ferentiable generalized synchronization 

The case of a non-differentiable synchronization map- 
ping formed in drive-response systems has been studied 
recently in In this paper, the stability of the re- 

sponse behavior in the whole phase space of the response 
system is linked with the existence of a mapping which 
maps the trajectories of the chaotic attractor in the driv- 
ing system into the trajectories of the response system 
with a mapping which is a Holder continuous function. 
These results help fill in the gap in parameter space which 
exists in between the synchronization regime where a 
smooth synchronization manifold exists and the natu- 
ral borders of the synchronization zones. In the present 
paper, we continue by studying how this theory of con- 
tinuous functions can be extended to cases of chaos syn- 
chrony with frequency ratio other than one-to-one, where 
the synchronization mapping can become multi-valued. 



In section ||, we present an example of chaos synchro- 
nization with frequency ratio 2:1 and study the forma- 
tion of the multi-valued mapping in detail. We do this 
by examining the bifurcations responsible for the transi- 
tion from the single-valued mapping to the multi-valued 
mapping as the regime of generalized synchronization of 



chaos changes. In section III, in order to provide a more 



complete view of the synchronization mappings and their 
properties, we supplement our numerical simulation re- 
sults with a discussion of some recent rigorous mathe- 
matical results not easily available to the reader, which 
are related to the considered regimes of generalized syn- 
chronization. We then review the results of the theory 
behind generalized synchronization with continuous func- 
tions and develop a similar theory for the case of multi- 
valued functions. 



II. CHAOS SYNCHRONIZATION WITH 
FREQUENCY RATIO 2:1 

We consider a regime of chaos synchronization with fre- 
quency ratio 2:1. This type of synchronization was ear- 
lier observed in the experiments with two directionally- 
coupled chaotic circuits. The sketch of the experimental 
setup is presented in Fig. |l|. The details on the exper- 
iment and the parameters of the circuits can be found 
in |0). 




Driving System 



Response System 



FIG. 1. The sketch of the circuit diagram for the exper- 
imental study of chaos synchronization with frequency ratio 
2:1. The frequency ratio is controlled by the multiplier used 
to select the values of the inductors and capacitors in the drive 
and response circuits. 

In this paper we present a detailed analysis of this syn- 
chronization regime based on numerical simulation re- 
sults. 



A. Model 

The dynamics of the drive circuit are described by the 
set of differential equations of the form El 
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VX2 — —Xi — 5X2 + 

vxz = 7(ai/(a;i) - xs) - ax2- 
The response system equations are 



(1) 



yi = 2/2 

y'2 = -yi - 5y2 
ys = i{a2f{yi) 



ys 

ys+gxi) - ay2, 



(2) 



where 7 = VLC T/R2C2 = 0.294, a = C1/C2 = 11.52, 
6 = m^/dT/L = 0.534, ai = 15.93 and = 16.7 
are the fixed system parameters and g is the couphng 
strength. Variables xi, X3, yi, and ?/3 correspond to the 
voltages across the capacitors (see Fig. |l|). Variables X2 
and 2/2 are proportional to the current through the induc- 
tors in the drive and response circuits respectively. The 
nonlinear function /( ) models the input-output charac- 
teristics of a nonlinear converter (N) used in the circuit. 
The shape of the nonlinearity is presented in Fig. ^ In 
the numerical simulation, we modeled the function f{x) 
with the formulae 



f{x) = sign(a;) (a~ ^ d{fp{x) - a)'^ + c\ /d 



where 



fp{x) 




if |a;| < a 
-p) if a < |x| < 6 
if \x\ > b 



d = (a"^ - c)/a^, q = 2a/{b - a), p = (6 + a)/2. The 
values of the parameters a, b and c are chosen to be equal 
to 0.5, 1.8, and 0.03 respectively (see |2^] for details). 
The parameter v in the equations of the drive system is 
the time scaling parameter, which is used to select the 
frequency ratio of the synchronization zone. 








FIG. 2. The shape of the nonlinear function /( ) measured 
in the experimental setup. 



The chaotic attractors occurring in the uncoupled drive 
and response systems with these parameter settings are 
shown in Fig. ^ Since we are interested in studying the 
chaos synchronization regime with a frequency ratio of 
2:1, the parameter u was set equal to 0.498. As a result, 
the phase velocity of the trajectories of the driving at- 
tractor is about as twice as high as that of the response 
attractor. 



(a) 



(b) 



yi 



FIG. 3. Chaotic attractors of the drive system computed 
with the parameter value ai = 15.93 (a), and the uncoupled 
response system (b). The bold dots on the trajectories of the 
driving attractor show the points where the chaotic trajecto- 
ries cross the Poincare cross section {x2 = with dx2/dt > 0). 



B. Onset of generalized synchronization 

Introduction of sufficiently strong coupling between 
the systems results in the onset of synchronization. The 
regime of synchronization which occurs when g = 2.5 is 
presented in Fig. ^. The onset of synchronization in this 
case is detected using Lissajous figures, by observing the 
attractors in the Poincare cross section, and by utilizing 
the auxiliary system method. Each of these methods will 
be briefly explained in the following paragraphs. 

The Lissajous figure is a standard approach frequently 
utilized in traditional studies of synchronization between 
periodic oscillators. It enables one to clearly see the onset 
of phase locking and define the frequency ratio of the 
particular synchronization zone. For the frequency ratio 
2:1, the Lissajous figure has a form similar to the shape of 
the digit "eight" . Despite the fact that due to the chaotic 
behavior the Lissajous figure does not appear as a closed 
curve, it is still not difficult to see the onset of phase 
locking and identify the frequency ratio (see Fig. ^). 

Analysis of generalized synchronization of chaos with 
the help of Poincare cross sections was proposed in . 
In this analysis, we define the Poincare cross section in 
the drive system and detect the moments of time when 
driving trajectory crosses it. At these moments, we sam- 
ple the state of the response system and compare the 
points to the corresponding set of the points in the drive 
system. This method simplifies the analysis of chaotic 
attractors by reducing their dimension by one. It helps 
to detect the onset of phase locking and, in some cases. 
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to examine the onset of topological equivalence between 
the attractors in the synchronized drive and response 
systems. In the numerical simulations we analyze the 
properties of the chaotic attractors on the Poincare cross 
section {x2 = 0) by computing the points where the tra- 
jectories cross the value X2 = with positive values of 
dx2 /dt (see Fig 0a) . 
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FIG. 4. (a) The Lissajous figure for the synchronized 
chaotic oscillations plotted on the plane of variables from the 
drive and response {xi,yi), and (b) the projections of the 
corresponding chaotic trajectories onto the plane of similar 
variables in the response and auxiliary systems {yi,zi). The 
strength of the coupling is g — 2.5. The bold dots show the 
points on the chaotic trajectories corresponding to the mo- 
ments of time when the trajectory of the drive system crosses 
the Poincare cross section {x2 — 0) with a positive value of 
dx2/dt . 

The auxiliary system method was proposed as a prac- 
tical method of detecting the formation of a continuous 
mapping through the ability of persistently point out the 
current state of the response system without direct com- 
putation of the map |Q. This method assumes the use 
of a third (auxiliary) system which is an exact replica of 
the response system. In our case, the auxiliary system is 
of the form 



Zl = Z2 

Z2 = -Zl - SZ2 + Z3 

Z3 = 7(02/(2:1) - 23 + 52:1) 



(3) 



aZ2- 



Due to the identity of the dynamics of the response and 
auxiliary systems, the nine dimensional space encompass- 
ing all three coupled systems has an invariant manifold 
y — z, where y and z are vectors of the variables in the 
response and auxiliary systems, respectively. When this 
manifold contains the entire chaotic attractor, a continu- 
ous mapping exists which projects the trajectories of the 
drive system onto the trajectories of the response systems 
after transients die out (see |^ for details). Therefore, 
the formation of a chaotic attractor in the invariant man- 
ifold (y — z) constitutes the onset of generalized synchro- 
nization between the drive and response systems. Such 
a chaotic attractor occurs in our simulation with g ~ 2.5 
(see Fig |b). 



When one utilizes this method in numerical simula- 
tions, it can be difficult to ensure that the chaotic set 
of trajectories in the invariant manifold does not contain 
transversely unstable limit sets. Due to the identity of 
the response and auxiliary systems and the limited accu- 
racy of numerical simulations, the chaotic trajectory can 
settle down on the invariant manifold y = z, even if the 
chaotic set in the manifold contains transversely unsta- 
ble limiting subsets. This artifact of numerical simulation 
can be resolved by breaking the symmetry between the 
response and auxiliary systems. In our simulations, we 
accomplish this by using slightly different values of cou- 
pling. The coupling of the auxiliary system is larger than 
the value of coupling in the response system by a value 
of Ag = 0.01, less than 1% of the coupHng strength used. 



C. Two regimes of generalized synchronization 

The goal of this section is to examine in detail the 
main features of two qualitatively distinct regimes of 
generalized synchronization observed for different val- 
ues of coupling in some previous experimental studies 
of chaos synchronization with a frequency ratio of 2:1 
(see [|o| for details on the experimental setup). In order 
to point out the change between the regimes of gener- 
alized synchronization versus variation of the coupling 
parameter, we analyzed the g dependence of dmaxis) = 

{\J (z — y)"^(z — y)}max, which is the maximal deviation 
between the response and auxiliary systems computed for 
the chaotic trajectories after transients die out. This de- 
pendence is presented in Fig. |^. Synchronization takes 
place for values of g where d„iaxig) is approximately 
equal to zero. In Fig. |^ there are two intervals of g which 
correspond to two synchronization regimes. These inter- 
vals are separated by a region of asynchronous behavior 
appearing around g — 2.0. 

The oscillations occurring in the synchronization inter- 
val g > 2.1 correspond to the chaotic attractor shown in 
Fig. ^. When the value of coupling decreases and arrives 
at the border of this synchronization interval, general- 
ized synchronization terminates. The chaotic attractor 
corresponding to these asynchronous oscillations, com- 
puted for g = 2.0, is presented in Fig. |^. The destruction 
of the generalized synchronization follows from the ap- 
pearance of sporadic outbursts of non-identical behavior 
in the response and auxiliary systems, which are seen in 
the Fig. ||b. Such outbursts are indicative of the exis- 
tence of transversely unstable limiting trajectories in the 
chaotic set located in the manifold y = z and therefore, 
indicates the inclusion of conditionally unstable limiting 
trajectories in the chaotic attractor residing in the phase 
space of the drive and response systems. 
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FIG. 5. The maximal deviation between the states of the 
response and auxihary systems versus the value of coupling 
parameter g computed with Ag = 0.01. 
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FIG. 6. The regime of asynchronous chaotic oscillations 
computed for g — 2.0. (a) The Lissajous figure in the plane 
{xi,yi), and (b) the projections of the corresponding chaotic 
trajectories onto the plane {yi,zi). The bold dots show the 
points on the chaotic trajectories corresponding to the mo- 
ments of time when the trajectory of the drive system crosses 
the Poincare cross section {x2 — 0) with a positive value of 
dx2/dt. 

Further decrease of the couphng strength again results 
in the onset of synchronization when the value of cou- 
phng is in the second interval of synchronization which 
takes place between g « 1.5 and g « 1.9 (see Fig. |^). The 
synchronized chaotic attractor which occurs in this case 
is shown in Fig. |^. In the joint phase space of the drive 
and response systems this attractor represents a closed 
ribbon formed after a "period-doubling" of the synchro- 
nized chaotic attractor studied in the first interval of syn- 
chronization, compare Figs. ^ and 0a. Due to the pe- 
riod doubling which occurs in the response system (the 
driving system remains unchanged), the synchronization 
mapping for this attractor becomes a one-to-two map- 
ping. Depending on initial conditions, the synchronized 
response oscillations settle down on the ribbon with one 
of two different phases. As a result, the single synchro- 
nized attractor in the drive-response system is presented 
as two different attractors in the drive-response-auxiliary 
system, one in the synchronization manifold y — z and 
the other one outside the manifold (see Fig. ^). 
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FIG. 7. The regime of synchronized chaotic oscillations 
computed for g = 1.8. (a) The Lissajous figure of the syn- 
chronized chaotic attractor in the plane of variables from the 
drive and response systems (a; i, yi). (b) Two chaotic attrac- 
tors corresponding to the synchronized chaos formed in the 
response and auxiliary systems plotted in the plane 
The bold dots show the points on the chaotic trajectories cor- 
responding to the moments of time when the trajectory of the 
drive system crosses the Poincare cross section {x2 — 0) with 
a positive value of dx2/dt. 

In order to examine the synchronous chaotic attractors 
in more detail and understand the properties of the cor- 
responding synchronization mappings, we study the lo- 
cation and stability of the limiting sets contained in the 
attractors. These limiting sets are the unstable periodic 
orbits (UPOs) embedded in the chaotic attractor. 



D. UPOs and the synchronization mappings 

It is well known that phase locking of the unstable 
periodic orbits in the coupled chaotic systems plays a 
significant role in the understanding of different regimes 
of chaos synchronization (see, for example fl^ , !!^ , ^ -^). 
Any UPO is characterized by the spectrum of multipli- 
ers or the corresponding Lyapunov exponents A^ = 
\n{^i)/Tupo, where Tjjpo is the period of the UPO. In 
the analysis of synchronization in directionally-coupled 
systems, it is always useful to split the whole spectrum 
of Lyapunov exponents into two groups: the first is given 
by the dynamics of the autonomous drive system, and 
the second comes from the conditional dynamics of the 
response system. The conditional dynamics corresponds 
the behavior induced in the response system by the driv- 
ing signal. 

Since the drive system is a three-dimensional dissipa- 
tive system in the present case, the first group of the Lya- 
punov exponent spectrum for each UPO contains three 
exponents: Ai > 0, A2 = 0, and A3 < 0. These expo- 
nents do not depend upon the coupling strength or any 
of the other parameters of the response system. The sec- 
ond group of exponents in the spectrum is due to the 
dynamics of the response system. Following Q we call 
these exponents the conditional Lyapunov exponents Af. 
If all of the conditional Lyapunov exponents of a UPO 
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have negative values, we call that orbit a Stable Response 
(SR) UPO. If at least one conditional Lyapunov exponent 
of the UPO is positive, we call the orbit an Unstable Re- 
sponse (UR) UPO. 
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FIG. 8. The Poincare cross sections of stable response (SR) 
UPOs, unstable response (UR) UPOs and chaotic trajectories 
in the drive-response system computed for (a) g — 2.5, (b) 
g = 2.0, and (c) g = 1.8. 



In our study, we consider all UPOs in the chaotic at- 
tractor of the driving system (shown in Fig. ^) up to 
period six, one of the period-8 UPOs, and one of the 
period-10 UPOs. Using the waveform of xi{t) of these 
orbits as a driving force, we computed the corresponding 
orbits that are formed in the driven response system. It 
is clear that, in this case, the stable periodic orbits corre- 
spond to the SR UPOs and the unstable ones to the UR 
UPOs formed in the joint phase space of the drive and 
response systems. 

The results are shown in Fig. ^, where we plot the 
intersections of the UPOs with the Poincare cross section 
{x2 = 0, dx2/dt > 0) projected onto the variable plane 
{xi, yi). In order to see if the UPOs belong to the chaotic 
attractor we also plot, in the background, the points on 
the Poincare cross section of the chaotic trajectory. 

When the systems are synchronized with g = 2.5, all 
the SR UPOs are inside the chaotic attractor, while the 
UR UPOs do not belong to the attractor, see Fig. The 
trajectory of the synchronized chaotic attractor wanders 
around the SR UPOs which form the skeleton of the at- 
tractor. Since this skeleton does not contain UR UPOs, 
the response system always follows the driving chaotic 
trajectory in a stable manner. This fact explains the sta- 
bility of the identical chaotic oscillations in the response 
and auxiliary systems (see Fig |^b) . 

Decreasing the value of the coupling parameter g re- 
sults in a sequence of bifurcations associated with the 
SR UPOs. These bifurcations create a large number of 
new UPOs, including new UR UPOs, which are added to 
the skeleton of the chaotic attractor (see Fig. ||b for the 
case g = 2.0). Since these new UR UPOs are inside the 
chaotic attractor they form regions where the response 
behavior for the chaotic driving is unstable. As a result, 
the synchronization is destroyed, see Fig. 

There is the coexistence of two types of UPOs which 
have different numbers of unstable directions in the 
chaotic attractor shown in Fig. |8|d. This is indicative 
of the onset of a non- hyperbolic situation [^7|-^9|. Note 
that even though some kind of non-hyperbolic situation 
may exist in the chaotic attractor of the driving system, 
it does not make any impact on the response behavior 
and therefore, makes no impact on the synchronization 
between the drive and response systems. The synchro- 
nization is sensitive only to the non-hyperbolicity of the 
chaotic attractors contributed by the response system dy- 
namics. This kind of non-hyperbolic situation is respon- 
sible for the inability to compute the trajectory of the 
response system if arbitrarily small perturbations of the 
system are taken into account This is the reason 

why the identical chaotic oscillations in the response and 
auxiliary systems do not follow the same path when even 
an arbitrary small perturbation are considered. It also 
directly indicates the destruction of the synchronization 
mapping. 
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FIG. 9. The bifurcation diagram yi{t„) vs g computed for 
the response system driven by the chaos (a), and by the 
UPOs: p2-{h), p4-(c), p6-(d), and p8-(e), which are found 
in the chaotic attractor of the driving system. 

As one can see from Figs. ^ and 0a, the further de- 
crease of the coupling between the systems results in the 
period doubling of the ribbon formed by the chaotic at- 
tractor. When this period doubling "bifurcation" is com- 
plete, all the SR UPOs and the chaotic trajectories of the 
attractor are moved far enough from the region occupied 
by the UR UPOs and, as a result, all UR UPOs again re- 
side only outside the chaotic attractor (see Fig||c). This 
qualitative change of the chaotic behavior eliminates the 
non-hyperbolic situation caused by the dynamics of the 
response system and, as a result, the generalized synchro- 
nization is regained (see 0b) . It follows from the analysis 
of the SR UPOs that in this regime of chaos synchro- 
nization any point on the UPOs contained in the chaotic 
attractor of the driving system maps in to two points 
located on the corresponding SR UPOs. Therefore, this 
synchronization is characterized by one-to-two synchro- 
nization mapping. 

In order to understand the properties of this synchro- 
nization mapping we studied the bifurcations of the SR 
UPOs when the value of coupling parameter is changed 
from high to low in the interval where the ribbon of the 



chaotic attractor experiences the period doubling bifur- 
cation. This study, plotted in Fig ||, reveals a very inter- 
esting fact. It turns out that not all the UPOs experience 
a period-doubling bifurcation in conjunction with the bi- 
furcation of the ribbon. Indeed, the orbits which have 
a period which is a multiple of four (the period of the 
chaotic ribbon after the bifurcation), are not subject to 
the period-doubhng bifurcations (see Fig ^ and e). The 
other UPOs go through the period doubhng bifurcation 
along with the ribbon. 
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FIG. 10. The bifurcation diagram yi{tn) vs g computed for 
the response system driven by the unstable periodic orbit p4 
located in the chaotic attractor of the driving system. The 
stable response orbits are indicated by circles and diamonds. 
The unstable response orbit is indicated by crosses. 

The difference in the bifurcation scenarios between the 
p(4xN) UPOs (with N an integer) and the other UPOs 
result in different synchronization mappings occurring 
when the response system is driven by these different 
UPOs. After the period doubling bifurcation, the or- 
bits whose periods are not multiples of 4 have a one-to- 
two synchronization mapping (see Fig. ^ and d), while 
the p(4xN) orbits do not bifurcate and their synchro- 
nization mapping remains one-to-one (see Fig. ^ and 
e). This situation raises the following question: How 
do these one-to-one synchronization mappings for the 
p(4xN) orbits fit into the one-to-two synchronization 
mapping which occurs for all trajectories of the chaotic 
attractor? The answer to this question comes from the 
analysis of the stable response images to periodic driv- 
ing with the p(4xN) orbits. When the ribbon bifurcates, 
the p(4xN) orbits respond with the creation of a new 
pair of response orbits. One is SR UPO and the other 
is UR UPO (this bifurcation scenario for the p4 orbit is 
presented in Fig. |lo|). Therefore, despite the fact that 
periodic synchronization for p(4xN) is characterized by 
a one-to-one mapping, each point on the p(4xN) orbit 
has two stable images in the phase space of the response 
system. Having two SR UPOs for each p(4xN) orbit, 
the orbits are able to conform to the one-to-two chaos 
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synchronization mapping. 

This example allows us to draw the following conclu- 
sions on the properties of generalized chaos synchroniza- 
tion with a multi- valued mapping. 

• The synchronized chaotic attractor represents a rib- 
bon of chaotic trajectories which does not con- 
tain limiting sets with unstable response dynamics. 
Each point in the chaotic attractor of the driving 
system maps into M points in the chaotic attractor 
in the joint phase space of the drive-response sys- 
tem (in the presented example, M — 2). The value 
of the number M is invariant over the attractor. 

• These M points are located in the chaotic ribbon in 
a such way that they represent the different phases 
of the motion along the ribbon. Any trajectory 
of the synchronized attractor passes through all 
phases in sequential order. Therefore, the multi- 
valued synchronization map can be presented as 
the cyclic sequence of the single- valued continuous 
maps. 

• The UPOs of the chaotic attractor in the drive sys- 
tem have SR UPOs, each of which is characterized 
by the mapping where one point on the driving 
UPO maps into M points of the corresponding SR 
UPO. The exception to this rule includes only the 
UPOs of the driving chaotic attractor with periods 
that are multiples of the period of the chaotic rib- 
bon. Each point on these special UPOs can map 
into a single point of the corresponding SR UPO. 
However, these UPOs should therefore have M sta- 
ble response UPOs. 

III. FUNCTIONS IN THE THEORY OF 
GENERALIZED SYNCHRONIZATION OF 
CHAOS 

In this section we formulate the theoretical framework 
for the synchronization mappings which occur in those 
regimes of chaos synchronization similar to the one con- 
sidered on the previous section. The main goal here is 
to link the properties of the synchronization mappings 
with the stability characteristics of the response system 
behavior. 

As it was shown in the previous section, the use of a 
Poincare cross section enables one to simplify the analysis 
of the synchronization mapping by studying the dynam- 
ics of the maps. In a rather general case, the maps occur- 
ring in the drive and response systems can be presented 
as maps with a skew product structure of the form: 

x„+i = r(x„) (4) 

y„+i = GA;(x„,y„), (5) 



where x and y are vectors consisting of the variables in 
the drive, response and auxiliary systems respectively, 
and fc is a parameter depending on the strength of the 
coupling. 

Since the first studies of chaos synchronization be- 
tween non-identical systems ||^, synchronization was in- 
terpreted as the system behaving in such a way that for 
any orbit (x„,y„), the coordinate y„ is just a function 
of the x„, once transients die out. The first rigorous re- 
sults on generalized synchronization were obtained for 
cases where there exists a smooth manifold belonging 
to the graph of a function y — /i(x) which is stable in 
the transversal direction and which contains the chaotic 
attractor corresponding to the synchronized oscillations 
||l^,|8|. This approach relies on the theory of normally 
hyperbolic invariant manifolds. 

In the most of the studies of different regimes of gener- 
alized chaos synchronization, it is observed that in some 
region of A:- values, a function h seems to exist, but it is 
not smooth. As the result, the graph of /i is a compli- 
cated geometrical object (see for example jl^,^). In this 
section we discuss recent rigorous results on the proper- 
ties of non-smooth functions h and present new results 
for the case of multi- valued mappings. 

A. Results and examples concerning single-valued 
functions h 

To make the paper self-contained, we first briefly dis- 
cuss the rigorous results on the existence and continuity 
of the synchronization function h which were obtained in 
a previous work [ p^ . The paper demonstrates that the 
function h is Holder-continuous if the contraction rate in 
the response system is small and Lipschitz-continuous if 
it is greater than a critical value. Assume that x G 3?™ 
and y e 3?" , and that is continuous and F is a homc- 
omorphism (i.e. F is continuous and invcrtiblc with an 
inverse F^^ which is continuous). The dynamics in the 
joint phase space of the drive and response systems is 
determined by a map 0fc : (x„,y„) i-^ (x„+i,y„+i). 
We also assume that in the joint phase space there 
exists a basin of attraction B = B^y S 3?™"'"", i.e., 
4>k{B) C Int{B) for any fc e 5, where S' is a region in 
fc-space in which the systems (^) and (^) exhibit master- 
slave chaos synchronization. Without loss of generality, 
we assume that B^y = B^ x By, i.e., B^y is a rectangle, 
where B^ (resp. By) is a ball in x-space (resp. y-space). 

Denote by Ak the maximal attractor in B^y, i.e., Ak = 
n^o'?^fe (^ajy)- Because of our assumptions, 

lim |y„ - y„| = 0, (6) 

n—>oo 

where (x„,y„) = (/)'j?(xo, yo), (x„,y„) = 0^(xo,yo,) and 
(xo,yo), (xo,yo) are arbitrary points in B^y. 
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To avoid non-essential technicalities it was supposed 
in that the limit in (p|) is achieved monotonically, i.e. 



|y«+l - Yri+ll < Wr. 



(7) 



for any k Cz S and any n > 0. 

Let Ak,x '■= ^x-A-k be the image of Ak under the nat- 
ural projection to 3?™. 

Theorem lA ||l9|] . (Existence) Under assumptions 
^ and (j^, there exists a function h : Ak,x 5^" such 
that Ak — graph(h). 

We emphasize that it is the inversibility of F that en- 
sures that /i is a single-valued function. Continuity plays 
no role at this stage. 

Theorem 2 A [|l9| . (Continuity) Under the assump- 
tions of Theorem lA, the function h is continuous. 

To obtain more detailed characteristics of this func- 
tional dependence, additional assumptions were made. 
It was assumed that 



|y«+i - yn+i| < ki\y. 



n J n\ 



(8) 



where fci < 1. Of course, the parameter fci is a function 
of k. For the sake of simplicity, we assume that k\ = k. 
Thus, 

|yn+i - y«-Hi| < fc|yn - Ynl, < fc < 1. (9) 
It follows that 



|Gfe(x,y)-Gfc(x,y)| <A;|y-y| 



(10) 



for any (x,y), (x, y) G Bxy. 

Assumption implies that |y„ — y„| goes to zero 
exponentially fast, and this fact allows one to prove that 
h is Holder continuous provided that the functions F and 
Gk have good smooth properties, or provided at the least 
that they are Lipschitz-continuous. 

Considering the forward and backward dynamics of the 
driving system we assume the following properties: 



and 



|F(x)-F(i)| <7+|x-i| 



|F-i(x)-F-i(x)|<7-|x-i| 



(11) 



(12) 



where 7_,7-|- > 1. When Lyapunov exponents do ex- 
ist (when the dynamics are differentiable with additional 
suitable conditions), the quantities log 7+, log7_ play the 
role of the forward and backward greatest Lyapunov ex- 
ponents respectively, and log k plays the role of the con- 
ditional Lyapunov exponent. 

It is assumed that function G^ (x, y) of the response 
system (^ is Lipschitz continuous with respect to x, i.e. 
for any (x,y), (x,y) G Bxy, 



|Gfe(x,y)-Gfc(x,y)| <??|x-x| 
where 77 > 0. 



(13) 



Theorem 3 A ||T^. (Holder property) Under assump- 
tions (|7^-(f7^ the function h is Holder continuous, i.e. 
for any < a < a^, and x,x Cz Ak,x one has: 



|/i(x) - /i(i)| < 2p|x-i|' 



where 



1 



1 - 



log(7+7-) ' 
loff k 



(14) 
(15) 



and p > pc, where pc is the solution of the equation: 



7^ iogl,7-[- 

p= rP 



logfc-loglR, 



T^+^xp % °' "' log(7+7-) 
7+ — A; V log k 



Theorem 4 A (Lipschitz property) Under the 

conditions of Theorem 3 A and provided that 



0<k< — , 

7- 

the function h is Lipschitz-continuous, i.e. 



|/i(x) - /i(i)| <L|x-i| 



(16) 



(17) 



where L > Lr 



The considered theorems give a pretty clear picture of 
the complexity of the synchronization mappings which 
usually occur in different regimes of generalized synchro- 
nization of chaos. It follows from these results that 
even in the case of differentiable generalized synchroniza- 
tion 17 1, the change of the parameters toward the border 



of the synchronization zone will gradually reduce the rate 
of contraction in the response system. At some critical 
value of the contraction rate (see Theorem 3 A and 4A), 
the smoothness of synchronization mapping will be de- 
stroyed while the systems remain synchronized with a 
non-differentiable function. 

Example: To confirm the validity of results above, 
let us consider coupled Henon maps. 

We consider the linear coupling of two non-identical 
Henon-type maps. For certain values of coupling 
strength, the hypothesis of our theorems are fulfilled. 

Let / be the following invertible map: / : [0,1]^ 
[0,1]2, /(xi,yi) = {x'^,y[) with 



Vi 



Vi 

hiVi) 



bxi 



where fi : [0, 1] — > [0, 1] is Lipschitz, with Lipschitz con- 
stant 7i, < 6 < 1. Let fi be such that f\{yi) < 
^ ,Vyi e [0,1], and (1 + 61)^ < 1. The map / is a 
homeomorphism of the unit square. Let gc be the foUow- 



mg map: gc 
with 



[0,1]'' [0,1]^ gcixi,yi,x2,y2) = (4,2/2) 
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X2 = y2 + c(yi - 7/2) 

y'2 = h{y2) + b2X2 + c((/i(2/i) + hixi) ~ (/2(2/2) + ^22:2)) 

where < c < 1, /2 : [0, 1] ^ [0, 1] is Lipschitz, with Lip- 
schitz constant 72, (l + fo2)-f- < 1 and /2(y2) < -.^Vi ^ 
[0, 1]. For convenience, we shall use the following short- 
hand: vi = and = (3:2,2/2 )■ 

To have contraction in W2-space, that is |gc(wi,i'2) — 
gciv\^V2)\ < k\v2 — W2I with < fc < 1, the following 
condition must be satisfied ^ 

1 



and 



c> 1- 



1 + &2 + f 



The forward and backward expansion rate of F are 
respectively given by: 7-|_ = 1 + ^ and 7^ —h^^. 

Theorem 3A holds, i.e. the function h is Holder con- 
tinuous and one can easily check the following expression 
for the Holder exponent, 



log[(l^c)(l + 62 + f] 

. / bl(l-c)(l + b2^) 

iog TT^I 



The condition ( [T6| ) in theorem 4 reads 

61 



c> 1 - 



1 + 62 + ? ' 



According to theorem 4A, this means that the function 
h is Lipschitz continuous. 



B. Multi- valued function h 

The framework of chaos synchronization discussed 
above applies only to the cases when the function h is 
a single- valued function, because of the assumption (|^). 
As a result this framework is not applicable to the case of 
synchronization presented in Fig. ^ where each point on 
the chaotic attractor Au.x in driving system maps into 
two different points of the synchronized chaotic attractor 
Ak- In this case By is not a simply connected region, and 
assumption (|^) cannot be used. Below we reformulate the 
results on single- valued function in order to extend it to 
the case of multi- valued function. 

Assume that By — Uf^^i?*, where are pairwise 
disjoint closed balls in the y-space, B^ is a closed ball in 
the x-space, and Bj^y = B^ x By. 

We also assume that monotonic synchronization occurs 



in B^ 



lim 



|yn 



(18) 



|yn+i -yn+i| < |yn -y«| 



(19) 



where (x„,y„) = (/'fc(xo,yo), (x„,y„) '/'fe(xo,yo,) and 
(xo,yo), (xo,yo) are arbitrary points in B such that yg 
and yo belong to the same ball for some i. 

Denote by Ak the maximal attractor in B^y^ i.e. Ak — 
^'^=n'K{'Bxy) and let Ak,x ■= ^xAk be the image of Ak 
under the natural projection H^; to 

In this situation some extensions of theorems 1-4 hold. 

Theorem IB . (Existence) Under conditions ( [7^ and 
^idj), and provided that for any x £ Ak,x and for any i, 
1 < i < p, there exist y^ g By such that (x, y^) G Ak, 
there exist p functions hi : Ak.x ~* i — ^, ■■■,P such 
that graph{hi) C B^ x By and Ak = U^^igraph(hi) . 

Scheme of the proof. The proof is mainly the same 
as the one for single- valued case (Theorem lA). How- 
ever we need to emphasize the following facts. For any 
X S Ak,x and for any n > 0, the set 0^(x, Li^^^By)) has p 
connected components inside (F"(x), Uf^^^i?* )) and every 
set (F"(x),i3*)) contains one and only one of them. In- 
deed, if we assume that (F"(x), By)) contains more than 
one connected components for some i then it will imply 
that there exist a j such that (F"(x),i3^)) contains no 
components. Since F is one-to-one, it means that there 
are no points of the attractor in (F"(x), B^)) and we have 
a contradiction. From this point forward, one can follow 
the proof of Theorem lA. □ 

Theorem 2B . (Continuity) Under assumptions of 
Theorem IB, the functions hi are continuous, for i — 
1, ...,p. 

Scheme of the proof. Let x and x be close to each 
other and consider the points (x, /i'(x)) and (x, /i*(x)). 
We show now that (j)'^^ (jx., k' (x)) and ^^^"'^(x, /i*(x)) be- 
long to the sets (F-1(x),S^) and (F^^Ix), B^)) , respec- 
tively, with the same number j. 

Assuming 

the contrary, i.e., (?!)^^(x, /i*(x)) G (F-1(x),S^) 
but (j)^\Sc,h\x)) G (F-i(i),sO) and / ^ j. 
Therefore 0fe(F^^(x), )) belongs to (x, ,5*) and 
0fe(F-i(x),BJ)) belongs to {x,,Bl) where i' ^ i. Con- 
sider now the points (F^^(x),y) and (F^^(x);y), y G 

Bl It follows that dzst((/.fe(F-i(x),y),(/.fe(F-i(x),y)^ 

is bounded from zero provided that x and x are close 
enough. Roughly speaking, this distance is greater 
than ^dist{B^y^B^y ). On the other side, we know that 
|F~^(x) — F~^(x)| is small, therefore by continuity of 0^, 
dist^(/)/j(F~-^(x), y) , ^/j(F~^(x), y)^ is also small, result- 
ing in a contradiction. 

The rest of the proof is the same as in Theorem 2A. 

□ 



*We use the norm \v\ — \x\ + \y\, where v = {x,y). 
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We assume now that conditions similar to (|8|)-([T3|) hold 
in Bxy, i.e., inequahties (|8|)-([To|) hold provided that y = 
yo, y = yo belong to the same connected component in 
the set By] other conditions remain the same. 

Taking into account the schemes of the proofs of The- 
orems IB and 2B, one can check that the proofs of the 
following two theorems are similar to the proofs of The- 
orems 3 A and 4 A. 

Theorem 3B . (Holder property) If the conditions are 
satisfied, under assumptions of Theorem IB, the func- 
tions hi, i = are Holder continuous, i.e. for any 
< a < ac, x,x G Ak.x one has: 



whe 



\h,{x) - h,{^)\ < 2p|x-x|' 



a < a 



' ^ _ log(7+7- 
loe k 



(20) 



(21) 



and p > Pc where pc is the solution of the equation: 

fj log(7_j_-y_) 



7+ 



■ P 



-+i 



xexp(A(fc,7+,7_)), (22) 



where 

A(fc,7+,7_) := 



logfc-log(max, 

r-H^ — log 7+7^ ■ 

logfc 



Theorem 4B . (Lipschitz property) Under conditions 
of Theorem SB and provided that 

1 



< fc < — 

7- 



the functions hi, i 
i.e. 



(24) 

p, are Lipschitz-continuous, 



|/i,(x)-/i,(x)| <i:|x-x| 



(25) 



where L > Lc '.— iTi-j ' 

In order to deal with problems in the case of multi- 
valued function h we make use of the auxiliary systems 
approach | |2^ . Consider an auxiliary system which is 
a replica of the response system (|^) and given by the 
equation of the form 



Zn+l — Gfc(x„, Z„) 



(26) 



Note that the auxiliary system (g6|) can serve, as well as 
in the case of standard generalized synchronization, to 
indicate the validity of synchronization. 

Let the dynamics in the joint phase space of the drive, 
response and auxiliary systems be determined by a map 
-ipk ■ (x„,y„,z„) i-> (x„+i,y„+i,z„+i). Due to the iden- 
tity between the response and auxiliary system there al- 
ways exist an invariant manifold Ms : y = z. Let us 



relate the ball to Bl: zq £ Bl iff yo = zq G B^. An 
indicator of multi- valued synchronization We can use the 
equality 



lim 



0, 



(27) 



where (x„,y„,z„) = ■(/'fc (xq, yo, zq) and (xo,yo,zo), is 
an arbitrary point in B such that yo S By and zq € Bl 
with the same index i, as an indicator of multi- valued 
synchronization. Indeed, it is simple to see that the fol- 
lowing statement holds. 

Theorem 5B . Under the assumptions of Theorem IB 
the equality (^^ is satisfied provided that yo G By and 
zo (z Bl with the same index i. 

The behavior of orbits in the attractor can be differ- 
ent. It could be related to a cyclic repetition, hi — > 
^2 ~* hi /i2 or a more complex sequence of 

hi. It depends on the partition of Ak,x into connected 
components. 

In the simplest case when Ak,x is connected, every or- 
bit in Ak behaves in the same way: for any i there is a j 
such that (l)k{graph{hi)) = graph{hj), i ^ i,j < p. This 
fact is a direct corollary of Theorem 2B. Indeed, for any 
two points Qi and Q2 in graph{hi) and any e > there 
is a collection of points Po, ...Pn G graph(hi) such that 
-Po = Qi, Pn Q2 and dist{P„P,+i) < e, < i < N-1. 
Because of the continuity of (j)k, the points (f>k{Pi) and 
(23) 4>k{Pi+i) belong to the same branch, say graph{hj) for 



any i — 0, ...N — 1. In fact, a trajectory switches the dis- 
joint balls in the y-space in a particular order, and this 
order is uniquely determined by a permutation i — > j. 

If Ak,x contains more than one connected components 
then the behavior of orbits inside each of them determines 
the same itinerary among branches graph(hi) . For exam- 
ple, if such a component contains the projection of a pe- 
riodic point then the itinerary will be also periodic. But 
if a component docs not contain a periodic orbit then the 
itinerary could be non-periodic; in this case Ak,x should 
have infinitely many connected components. 



IV. CONCLUSIONS 

The example of generalized synchronization of chaos 
considered in Section || enables us to explore properties 
of a new type of chaos synchronization mapping in which 
the mapping is a multi- valued function. Based on analy- 
sis of the conditional stability of the synchronous chaotic 
behavior of the response system, we come to conclusions 
concerning the existence and continuity of the synchro- 
nization mappings formed in these regimes of generalized 
synchronization. We used the auxiliary system method 
to detect the onset of conditional stability both in the 
case of single- valued mappings and in the case of multi- 
valued mappings. 
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Changing the strength of the couphng between the 
drive and response systems, we followed the transition 
from a regime of synchronization with single- valued func- 
tion to a regime with a double-valued function. These 
two regimes of synchronization are separated by a regime 
of asynchronous chaotic oscillations. Our numerical anal- 
ysis showed that this regime of asynchronous oscillations 
occupies an interval of the coupling parameter values 
where all of the necessary bifurcations occur which are re- 
quired for the formation of the new type of synchroniza- 
tion mapping. Analysis of the unstable periodic orbits 
(UPOs) contained in the chaotic attractor formed in the 
joint phase space of drive and response systems revealed 
a very interesting element of this bifurcation scenario. 
We have found that not all UPOs in the attractor ex- 
perience the period doubling bifurcations as the ribbon 
containing the chaotic trajectories of the attractor dou- 
bles its period. The UPOs with periods that are multiples 
to the period of the post-bifurcation chaotic ribbon re- 
main qualitatively unchanged. In this case the response 
system generates additional UPOs with identical periods 
via tangential bifurcations. Thanks to these additional 
UPOs the newly formed double-valued synchronization 
function applies universally to all orbits of the driving 
chaotic attractor. 

The introduction of multi- valued synchronization map- 
ping into the concept of the generalized synchronization 
of chaos provides a new theoretical framework which is 
crucial for understanding of chaos synchronization phe- 
nomenon in many cases of synchronization in direction- 
ally coupled systems. Such cases include synchroniza- 
tion with frequency ratio other then one-to-one where 
the formation of multi- valued synchronization mappings 
is a quite natural occurrence. 
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